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We use new kinematic data from the ultra-faint Milky Way satellite Segue 1 to model its dark 
matter distribution and derive upper limits on the dark matter annihilation cross-section. Using 
gamma-ray flux upper limits from the Fermi satellite and MAGIC, we determine cross-section ex- 
clusion regions for dark matter annihilation into a variety of different particles including charged 
leptons. We show that these exclusion regions are beginning to probe the regions of interest for a 
dark matter interpretation of the electron and positron fluxes from PAMELA, Fermi, and HESS, 
and that future observations of Segue 1 have strong prospects for testing such an interpretation. We 
additionally discuss prospects for detecting annihilation with neutrinos using the IceCube detector, 
finding that in an optimistic scenario a few neutrino events may be detected. Finally we use the 
kinematic data to model the Segue 1 dark matter velocity dispersion and constrain Sommerfeld 
enhanced models. 
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I. INTRODUCTION 

The Fermi Gamma-ray Space Telescope and the cur- 
rent generation of Imaging Atmospheric Cherenkov Tele- 
scopes (ACTs) such as VERITAS, MAGIC, HESS, and 
CANGAROO, are exploring the gamma-ray sky to un- 
precedented accuracy. Additionally, neutrino detectors 
such as IceCube will be probing the neutrino sky with 
greater sensitivity, lower energy thresholds, and finer an- 
gular resolution than achieved to date. If weakly interact- 
ing massive particles (WIMPs) constitute the dominant 
component of the dark matter in the Universe, experi- 
ments such as Fermi, ACTs, and IceCube may be sensi- 
tive to their annihilation or decay products. 

Dark matter dominated Milky Way (MW) satellite 
galaxies are excellent targets to search for dark mat- 
ter annihilation products, because they are nearby and 
largely free of astrophysically-produced high energy pho- 
tons. Among the population of MW satellites, of par- 
ticular interest is the ultra-low luminosity object Segue 
1, which was discovered in 2006 as an overdensity of re- 
solved stars in the Sloan Digital Sky Survey (SDSS) [T] 
at a heliocentric distance of 23 ± 2 kpc. Originally clas- 
sified as a tidally-disrupted globular cluster, follow-up 
spectroscopy for individual stars has revealed that Segue 
1 is the least luminous and most dark-matter dominated 
galaxy known 2J. Though its proximity to the Sagit- 
tarius stream has led to a suggestion that Segue 1 is a 
dissolved star cluster originally associated with the Sagit- 
tarius dwarf galaxy j3j, new measurements of Segue 1 
member stars strongly support that it is indeed a dark 
matter dominated object [11 [5]. 

Given its location at 23 ± 2 kpc and the fact that it 
is dark matter dominated. Segue 1 is an ideal source for 
dark matter studies [H |3]. Since Segue 1 can be eas- 



ily localized, any signal or limit on high energy gamma- 
ray/neutrino emission from it is more straightforward to 
interpret than a corresponding diffuse signal from the 
Milky Way halo. As we show in this paper, the integral 
over the square of the dark matter density distribution is 
similar to, and potentially larger than, the same quantity 
for all other known Milky Way satellites. 

In this work we give the cross-section exclusion regions 
for dark matter annihilation into a variety of Standard 
Model particles using gamma-ray flux upper limits for 
Segue 1. We also give exclusion regions for dark matter 
annihilation to charged leptons, where an accompanying 
gamma-ray signal is guaranteed to exist from "final-state 
radiation" (FSR) [SHE]- The resulting flux from FSR pho- 
tons is independent of the astrophysical environment in 
which the photons are produced (as opposed to inverse 
Compton or synchrotron processes). Neutrino produc- 
tion is also guaranteed from the decay of muon or tau 
leptons, and we also investigate the case in which dark 
matter annihilates purely into neutrinos. 

Constraints on dark matter annihilation to charged 
leptons are particularly relevant in light of the positron 
results from the PAMELA satellite in the energy range 
of 10-100 GeV [H] as well as the Fermi [TO] and HESS 
[in [12] e+ -h e" spectra at ~ 0.3 - 1 TeV. If dark matter 
annihilation is responsible for these measurements, the 
annihilation cross-section must be 0(100 — 1000) times 
larger than the thermal freeze-out cross-section and only 
a small number of anti-protons can be created in the an- 
nihilation process [13n20 . Given the importance of the 
dark matter interpretation of these measurements, and 
the fact that more precise e+/e^ flux measurements may 
not conclusively determine their origin (given the uncer- 
tainty in the e^/e~ propagation in the Galaxy and in 
their background fiuxes) , independent checks of the dark 



matter interpretation of the PAMELA, Fermi, and HESS 
signals are of the utmost importance. 

Here we present annihilation cross-section constraints 
from the 9-month Fermi gamma-ray flux upper limits for 
Segue 1. We also present constraints from gamma-ray 
flux upper limits obtained by MAGIC. If dark matter an- 
nihilation is responsible for the electron and/or positron 
signals observed by PAMELA and/or Fermi/HESS, we 
also show that the ACTs MAGIC and VERITAS have 
excellent prospects for testing models and detecting dark 
matter by observations of Segue 1. Finally, we investi- 
gate the prospects for detecting a neutrino signal from 
dark matter annihilation in Segue 1 with IceCube, and 
we give constraints on Sommerfeld enhanced models. Re- 
garding Sommerfeld models that explain both Fermi and 
PAMELA, we find that current limits slightly disfavor 
light mediators of mass 0(0.1 GeV) if we take the mean 
value for the Segue 1 line-of-sight integral of the dark 
matter halo density-squared. 



II. FLUX CALCULATION 

The gamma-ray/neutrino flux from annihilating dark 
matter is 
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where m^ is the dark matter particle mass, {av) is the 
annihilation cross-section, Eth is the threshold energy of 
a given gamma-ray/neutrino instrument, and E^ax is the 
maximum energy of the photons/neutrinos. The integral 



over 1^ depends only on the particle physics details of 
the dark matter annihilation. £ann is the line-of-sight 
integral over the square of the dark matter density 
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and depends only on the properties of the dark mat- 
ter halo and the solid angle over which it is observed. 
Here p{r) is the halo dark matter density profile, r = 
\/D'^ + s^ — 2sDcos6 for a halo at a distance D, and the 
integral is performed along the line of sight over a solid 
angle Afi = 27r(l — cos 6*). Here we assume (av) is con- 
stant throughout the halo. 

Dark Matter Distribution — To determine £ann we 
use the new sample of line-of-sight velocities of Segue 1 
member stars from [3J [S]. This sample consists of 71 
probable member stars (as compared to the sample of 
24 members used in previous indirect detection stud- 
ies [2| [2TH24] ). In our calculations we include 66 stars, 
excluding from the sample of 71 two RR Lyrae variable 
stars and two probable binaries identified in [1] , and also 
excluding the star SDSSJ100704.35-hl60459.4 which is a 
6ct outher from the mean velocity of Segue 1 [98]. (See [4] 
for a detailed discussion of these stars and the sample in 
general, including the effects of binary star contamination 



Within a radius of i 
(degree) 



A-ann 

logiJ£a„n/(GeV"cm-5) 



0.01 
0.1 

0.25 
1.0 



18.3 ±0.9 

19.0 ±0.6 

19.1 ±0.6 

19.2 ±0.6 



TABLE I: Segue 1 C^nn values within different angular re- 
gions, with 6 — 0.25° corresponding to the radius of the out- 
ermost confirmed member star. 



[5S].) We adopt the maximum likelihood procedure that 
is described in detail in [211 ES] , and review the aspects 
of the calculation that are important for the purposes of 
this paper. We assume a gaussian likelihood in the dis- 
tribution of stellar velocities, with a dispersion given by 
the quadrature sum of the measurement uncertainty of 
each stellar velocity and the intrinsic dispersion, the lat- 
ter of which is given by a solution to the Jeans equation. 
The Jeans equation takes as input the dark matter den- 
sity profile, stellar density profile, and the stellar velocity 
anisotropy. We assume that the dark matter and stel- 
lar density profiles are independent. For the dark matter 
density profile we assume a three-parameter Einasto pro- 



A profile of this 



file, pir) =p,exp{-f[(;^) - l] } 

form is motivated by numerical simulations, which sug- 
gest that a = 0.17 provides a good fit for dark matter 
halos over a wide variety of mass scales [27]. (Though 
we use the Einasto parameterization in this paper, the 
results we present would not change if we use a differ- 
ent parameterization, e.g. the generahzed NEW model 
considered in \'22\.) For simplicity in this paper we con- 
sider isotropic models with /3 — 0; we find that this re- 
striction has minimal impact on our resulting constraints. 
The Segue 1 three-dimensional stellar density profile used 
in the Jeans equation is taken as a Plummer profile, 
p{r) ex (1 ± r^/rp)"^/^, where Vp = 30 pc. We take 
uniform priors on the following parameters over the fol- 
lowing ranges: a € [0.14 : 0.3], Tg/kpc £ [0.01 : 0.6], and 
logio [Ps/(M0kpc^'^)] e [6 : 10]. We then marginalize 
over the gaussian likelihood to construct the likelihood 
for these parameters. Note that given our assumption of 
uniform priors the likelihood is equivalent to the poste- 
rior probability density distribution for the parameters. 
From the likelihood for the parameters ps, rs, and 
a, we construct the likehhood for £ann using Eq. Q. 
To remain conservative we integrate the density squared 
within a solid angle oi 9 = 0.25°, corresponding to 100 
pc. This is close to the radius of the outermost known 
bound star in the Segue 1 population, which lies at 87 
pc. We fit gaussians to the log of the £ann likelihood, 
and for Segue 1 we find 



C 



Scgucl 



10 



19.1±0.6 



GeV^cm~^ 
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with Icr error bars and 6 within 0.25°. See Table U for 



-Cann values of Segue 1 within other radii. 

For the calculations below of the expected Sommer- 
feld enhancement, we are also interested in the velocity 
dispersion of the dark matter. Segue 1 is a particularly 
interesting, if not the most interesting, galaxy to consider 
in the context of Sommerfeld enhanced models because 
of its extremely low velocity dispersion. In our calcu- 
lations we will consider the one-dimensional radial dark 
matter velocity dispersion, <Tr, averaged over the three- 
dimensional dark matter mass distribution. We define 
this quantity as 



J p{r)r'^dr 



(4) 



where p{r) is the Einasto profile as above. From a given 
density profile, the one-dimensional dispersion a^ is de- 
rived from the Jeans equation, assuming that the dark 
matter velocity dispersion is isotropic. The likelihood for 
this quantity may then be determined from our base set 
of three parameters in a manner similar to the derivation 
of £ann- We find 



VTCT|y^4.5±1.8km/s. 



(5) 



We stress again that this quantity is the radial velocity 
dispersion of the dark matter, which we find to have a 
mean that is slightly larger than the observed line-of-sight 
velocity of stars reported in Ref. 0] . 

The value of £ann"°^ can be compared with a canon- 
ical value from the Galactic center. Using an Einasto 
density profile with standard Milky Way fit parameters, 
Ps = 0.06 GeV/cm3 and r, = 20 kpc [Mj, we find 
/:MW ^ 1020-82 Qev2 cm-5 within a radius of 0.25°. Al- 
though this is about two orders of magnitude larger than 
the mean value of /Canf"''^ within the same solid angle, 
the large astrophysical background in the Galactic cen- 
ter makes Segue 1 a much cleaner target. Moreover, the 
Milky Way density profile is estimated from simulations 
that have not included the effect of baryons, whose im- 
pact on the density profile is not clear. In addition, 
these simulations assume that the dark matter is cold 
(e.g. a thermal WIMP or an axion), but a non-thermal 
WIMP, or a WIMP with non-negligible self-interactions, 
could give different halo profiles [29|. Though the Milky 
Way dark matter distribution can be probed by stellar 
kinematics on larger scales (e.g. [3D]), the distribution 
near the Galactic Genter (within less than ~ 100s of pc), 
where most of the annihilation flux emanates, is currently 
completely unknown. 

The mean value of £ann for Segue 1 is the same as for 
Draco |31j . and is larger than that for any other known 
dwarf. Although previous studies have derived larger 
'Cann valucs for Willman 1 and Sagittarius [55], neither 
system appears to be in dynamical equilibrium [351 [33] , 
which makes constraints on their dark matter halos from 
existing data unreliable. 

Gamma-rays from Dark Matter — For the direct 
annihilation channels into leptons, XX ~^ ^^i where i is 



a lepton, the energy spectrum from final state radiation 
(FSR) is given by 



dN-y _a fl + il-yy 
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[71I1II1I31H35]. Here a ~ 1/137, y ^ E^/m 



and Er, 



in Eq. M 



This formula holds in the 



collinear limit, where the photon is emitted coUinearly 
with one of the leptons and mi <C fn^ . 

For the channels that include r's, the photon contri- 
bution from T decays into tt'^'s (which then decay to 7's) 
dominates over the FSR signal, and the energy spectrum 
is given in [551 [SZ] ■ The energy spectrum for FSR for the 
channel xX — > <^</' — > (-M is given explicitly in [55]. If 
the final state includes /i's, we also include photons from 
the radiative decay of the muon, e.g. /z~ — )■ e~v^Ve'~i, us- 
ing the formulas found in e.g. [55] [351 ISH]- We do not 
consider gamma-ray signals from synchrotron radiation 
or inverse Compton scattering as both are expected to 
be negligible in Segue 1. 

For dark matter annihilation into W^W~ {W bosons), 
hb (bottom quarks), ti (top quarks), and gg (gluons), we 
use the photon spectra from DarkSUSY lU |42] . Pho- 
tons produced from these annihilation channels tend to 
be softer than FSR photons, but more numerous. This 
makes Fermi more sensitive than AGTs to these annihi- 
lation channels. For example, for the Wino LSP model 
[10] that explains the PAMELA e+ and p data (but 
not the Fermi or HESS e+ -t- e~ data), the dominant 
channel is xX ~^ W~^W~ , and the photon contribu- 
tion comes from W decays generating tt'^'s and hence 
7's. To compare with Fermi flux upper limits, the rel- 
evant integral over dNj/dE^ for this model is given by 
J^j^^y{dNy/dE^)dE^ ~ 27.14, where m^ ^ 180 GeV 
j40j . Since most of the photons have energies around 1 
GeV, ACTs will not be very constraining for this model 
even if they lower their energy threshold to ~ 50 GeV. 

Neutrinos from Dark Matter — We estimate the 
number of muons detected by IceCube from muon neu- 
trinos coming from dark matter annihilation in Segue 1 
following e.g. [321113] (see also [H]). We consider the di- 
rect channel xX ~^ ^/xJ^p) which gives a good indication 
of how well IceCube could do in a "best-case" scenario, 
as well as the two channels xX ~^ M^M" ^i^d xX ~^ 4'4> 
with <j) — > fJ.^ p^ (we do not consider channels with taus, 
although their decay would also produce neutrinos) . The 
muon decays p~ — >■ ^v^v^e~ and /i"*" — >■ ^v^v^e'^ produce 
both electron and muon (anti-)neutrinos. The neutrino 
energy spectra for the direct channel xX ^ ^'^i'/i is given 
by 
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dE, 



:^ = 2SiE,-m^), 



while for the channel xX ~^ A'^M we have 
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FIG. 1: Current exclusion regions from Fermi 9-month gamma-ray observations of Segue 1 (bounded below by tiie blue solid 
diagonal line) and MAGIC gamma-ray observations of Segue 1 (bounded below by the purple solid curved line). The exclusion 
regions use the conservative 2ct lower limit of Xann given in Table ^ within 6 = 0.25 (Fermi) and = 0.1 (MAGIC). The 
dashed blue and purple lines depict the respective cross-section bounds using the optimistic 2cr upper limits of Xann- For the 
XX ~^ W"^W~ channel, the black dot is the region favored by a model of wino-like neutralinos that explains the PAMELA 



positron data HOj. Note that 



mt 



175 GeV for the XX ~^ t^ channel. 



The energy spectra for the longer channel xX 



fj, II are given by 
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where x — E^/m~f^ (39] . (The Vi and Vi spectra are the 
same for each channel). When these neutrinos reach 
Earth, the probability that Vi will have oscillated into 
v^j_ is roughly |45| 
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P{i>e 






'. 0.39, 
0.22 



(12) 



(we take this oscillation probability also into account for 
the direct channel xX ~^ t/^^^)- Since Segue 1 lies in the 



Northern Hemisphere, these neutrinos travel through the 
Earth towards IceCube. While i>e and i^r predominantly 
give rise to cascade-like events in IceCube, I'^'s can con- 
vert to muons in the ice and produce track-like events 
that yield much better angular resolution [15]. We thus 
focus exclusively on detecting the muons from the i^^^'s. 
(In particular, we also ignore the muons produced from 

Vr ^ T ^ fl.) 

Given dN^ ^p /dEi, above, one obtains the differential 
neutrino flux, d^^^,jj^JdE^, from Eq. (|l|). The muon 
energy spectrum detected by IceCube in a time T is given 
by 
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FIG. 2: Current exclusion regions from Fermi 9-month gamma-ray observations of Segue 1 (bounded below by the blue solid 
diagonal hne) and MAGIC gamma-ray observations of Segue 1 (bounded below by the purple solid curved line). The exclusion 
regions use the conservative 2a lower limit of £ann given in Table^within = 0.25 (Fermi) and = 0.1 (MAGIC). The dashed 
blue and purple lines depict the respective cross-section bounds using the optimistic 2a upper limits of £ann- Green and orange 
regions indicate the best-fit regions to explain the PAMELA and Fermi e"'"/e~ signals and are adapted from |^ after rescaling 
from a local dark matter density of 0.3 GeV/cm^ to 0.43 GeV/cm^ [47] . 



Here the differential cross-section for a neutrino of energy E^ to scatter off a nucleon with mass ttin is roughly given 
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FIG. 3: Left: Projected number of muon signal events coming from neutrinos originating from dark matter annihilation 
directly into i^^i^,, in Segue 1, after 10 years of observation with IceCube, and assuming the mean value of £ann"'^^ given 
in Eq. (|3|. Right: Prospective exclusion regions for IceCube 10-year observations of Segue 1, assuming (signal muon 
events) /v signal + background muon events > 2. The muons observed by IceCube are assumed to come from muon neutri- 
nos produced in dark matter annihilation. Here the solid (dashed) line uses the 2a lower (upper) limit on £ann in Eq. S. 
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while for an anti-neutrino the numbers 0.20 and 0.05 are 
interchanged. For this expression, we took the average 
of the scattering cross-sections off a neutron and off a 
proton, and we assumed E^ ^ m^^/mN ~ 7 TeV, where 
mw is the W boson mass. (For dark matter masses of 
~ 10 TeV that give rise to neutrinos of ~ 10 TeV, the 
constraints we obtain will thus be optimistic by a factor 
of a few.) Once produced at an energy Efj,, a muon will 
loose energy according to 



dr 
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where a ~ 2x10"^ GeV cm^ g-i, /? ~ 4.2x 10"^ cm^ g-i 
[45] . and pm is the Earth matter density. If we require 
that the muon be detected above an energy threshold of 
Eth, it will travel a distance 
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The muon must intersect IceCube's effective area, which 
for E^ ~10 TeV is roughly 1 km^, but is much less 
for a lower energy muon. We parameterize the effective 
area, Aes{E^), for the full IceCube detector following 
[so t I51 j . To model IceCube's energy resolution, we also 
smear the muon energy spectrum given in Eq. ( 13 1 with 
logio(i?nieasurcd/£'roai) givcu by a gaussiau distribution 
centered on zero and with a standard deviation of 0.15 
[46] . The total number of muons detected by IceCube is 
then given by integrating this smeared energy spectrum 
over all muon energies from i^th up to m^. 



To calculate the muon background, we find 
dN^^^jj^/dEi, from atmospheric neutrinos using an 
average over zenith angle from |52j . We assume a con- 
stant pointing resolution of 1° radius down to neutrino 
energies of 100 GeV for IceCube ^53J, which should be 
obtainable given Segue I's zenith angle of ~ 106° with 
respect to the South Pole. 

By fixing the ratio of signal muons to the square root 
of the signal plus background, {av) can be determined as 
a function of m~^. 

While the experimental properties quoted above refer 
specifically to IceCube, the inclusion of DeepCore [46] in 
conjunction with IceCube will help in the reconstruction 
of muon events (ensuring a better pointing accuracy) and 
in the differentiation of background events. 



III. RESULTS 
Current and Projected Observations 

Below we discuss current and projected experimental 
flux limits and resulting cross-section bounds. 

Fermi — Fermi has reported gamma-ray constraints 
on 14 nearby dwarf galaxies from 11 months of data in 
[31) . For Segue 1, Fermi has only reported a flux limit of 
4x 10~^° cm~^ s~^ for photons above 100 MeV and below 
300 GeV using 3 months of data [MllSSj . They assumed a 



source spectrum of 



dE-, 



1 



which is very close to the 



FSR spectrum. (Note that the photons from r decays 
into 7r°'s follow a different spectrum.) Since no signal 
was seen from Segue 1 in the 9- month data [Ml |57] , we 
conservatively take the flux limit from the 3-month data 
and divide by v3 to arrive at an approximate 9-month 
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FIG. 4: Projected number of muon signal events coming from neutrinos originating from dark matter annihilation directly into 
muons (left) or into muons through an intermediate particle (p (right) in Segue 1, after 10 years of observation with IceCube, 
and assuming the mean value of /ifnl"'^^ given in Eq. (|3|. Green and orange regions indicate the best- fit regions to explain the 
PAMELA and Fermi e'^ /e~ signals and are adapted from [20 after rescaling from a local dark matter density of 0.3 GeV/cm^ 
to 0.43 GeV/cm^ [47] • Note that these projections do not assume any additional Sommerfeld enhancement in the dwarf galaxy 
compared to locally. However, if present, such an enhancement could effectively shift up by an order of magnitude or more the 
PAMELA/Fermi regions in the right plot, increasing the expected number of signal events. 



flux limit (2.3 x 10^^ cm^ ^ )■ (Fermi has not yet pre- 
sented flux limits on a l/E^ spectrum using more than 
3 months of data. The analysis in [23] uses 9 months of 
data to constrain the "Constrained Minimal Supersym- 
metric Standard Model".) This simple scaling expected 
from statistics is corroborated by Fermi's results for a 
1/-E^ spectrum for which they did present results us- 
ing 3 and 9 months of data: the flux limits for Segue 
1 were 3.1 x 10"^ cm^^ s'^ [SS] and 1.83 x IQ-^ cm^^ 
s~^ [5S1[S7|, respectively. (It does, however, ignore any 
improvements that may have been made to the analy- 
sis in the meantime, which could further strengthen the 
limits.) As expected, these limits are worse than the flux 
limit for a 1/ E^ signal spectrum, which has more photons 
at higher energies. 



ACTs — Segue 1 has been observed by MAGIC and 
a flux upper limit of roughly lO^^^cm"^ s~^ above 100 
GeV has been reported [SS] . This limit should be viewed 
as an order of magnitude estimate as the exact limit will 
depend on the shape of the spectrum considered. This 
limit is for observations of Segue 1 within 0.1 degrees of 
its center. 



IceCube — For IceCube, we show both the 
number of signal muon events expected and 
the cross-section bounds assuming (signal muon 
events)/-\/signal - 
assume an 



background muon events > 2. We 
energy threshold of 100 GeV, a pointing 
resolution of 1° in radius that is constant with energy, 
and 10 years of observations. 



Cross-Section Bounds 

Prom Gamma-rays — In Figure [ij we show the cur- 
rent upper bounds on the dark matter cross-section as 
a function of mass from Fermi and MAGIC observations 
of Segue 1. We consider the channels XX ~^ W^W~, 
bb, tt, and gg. Figure^ shows the current upper bounds 
for the channels xX ^^ e+e^, xX ^^ A^^M ? XX ^^ t^t, 
XX^ H ^ e+e"e+e", xX ^ # ^ M^A*"A*^M", and 
XX ~^ 4'4>-, with (j) — > e+e~ 90% of the time, and — > 77 
10% of the time. Here x denotes the dark matter par- 
ticle, while (j) denotes an intermediate particle such as a 
gauge boson or a scalar. For the case that (j) goes to 77 
10% of the time, is a light scalar that mixes with the 
Standard Model Higgs boson. The channels in Figure [2] 
can give rise to the e+ and e+ -I- e~ excesses observed by 
PAMELA and Fermi/HESS, as can various combinations 
of these channels with different branching ratios. In all 
cases, we consider the mass of 4> to be m^ < ©(GeV). 
This kinematically constrains the (p to decay to Stan- 
dard Model particles that include a sizeable lepton frac- 
tion and avoids the overproduction of anti-protons. This 
can also lead to Sommerfeld enhanced annihilation (see 
below). 

In these flgures, the solid lines correspond to con- 
straints using the mean £ann values minus the 2a er- 
rors as given in Table ll] within 6 — 0.25 (for Fermi) 
and 6* = 0.1 (for MAGIC). Thus the sohd lines show 
very conservative 95% confldence level upper limits. The 
dashed lines use the mean £ann values plus the 2cr er- 
rors showing an optimistic reach. The green and or- 
ange regions indicate the best-fit regions to explain the 
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FIG. 5: Prospective exclusion regions for IceCube 10-year observations of Segue 1. The IceCube constraints assume (signal 
muon events) /-^signal + background muon events > 2. The muons observed by IceCube are assumed to come from neutrinos 
produced by dark matter annihilation directly into muons (left) or into muons through an intermediate particle (p (right). The 
solid (dashed) line uses the 2a lower (upper) limit on £ann in Eq. ([3|. Green and orange regions indicate the best-fit regions 
to explain the PAMELA and Fermi e"'"/e~ signals and are adapted from [20] after rescaling from a local dark matter density 
of 0.3 GeV/cm^* to 0.43 GeV/cm^' pff. 



PAMELA and Fermi e+/e^ signals and are adapted from 
[20] . Note, however, that we scaled down the regions in 
[5D] from a local dark matter density of 0.3 GeV/cm'^ to 
0.43 GeV/cm3 ;47]. 

Figures [T] and [2] show that Fermi has an advantage 
over ACTs for those dark matter channels that pro- 
duce softer photons. Fermi's much lower photon energy 
threshold compared to that of the ACTs (~100 MeV ver- 
sus ~100 GeV) also gives it an advantage in determining 
cross-section bounds for lower dark matter masses 0(100 
GeV). ACTs have an advantage over Fermi in constrain- 
ing models which produce harder photon spectra, assum- 
ing that the dark matter mass is above the ACT's thresh- 
old of around 100 GeV. This is because ACTs in general 
have larger collecting areas than Fermi (^ 10"* m^ versus 
^ 1 m^) and thus have better statistics at high photon 
energies. 

In Figure [I] the results for the XX "^ W^W~ chan- 
nel show that the Wino LSP model can be probed by 
gamma-ray observations of Segue 1. This model, which 
can also explain the PAMELA data, predicts an annihila- 
tion cross-section of ~ 2.5 x 10"^"' cm^ s~^ for m-^ ~ 180 
GeV (indicated by the black dot) [ID] . While the obser- 
vations disfavor this model (requiring a value for £ann 
less than its mean for consistency) they are unable to 
rule it out due to the uncertainty in £ann- (Note that 
this model also predicts gamma-ray lines, and limits on 
such lines disfavor it [59]). A more precise determina- 
tion by the Fermi collaboration of the flux limit, using 
the photon spectrum for dark matter annihilation into 
W^W~ , may set a stricter limit. 

In Figure [2] we see that ACTs give stricter con- 
straints than Fermi on dark matter models favored by 
the PAMELA and/or Fermi/HESS cosmic-ray electron 



and positron data. In many cases, the mean value for 
'Cann probcs the relevant cross-section and mass region, 
although taking the conservative value for £ann only con- 
strains the XX ~^ T'^T~ channel. Below we discuss 
in detail how Sommerfeld enhanced annihilation can be 
probed with Segue 1, which can lead to stronger con- 
straints in the case of dark matter annihilation to (/)'s. 

Prom Neutrinos — In Figure [s] (left) , we show the 
number of signal muon events expected as a function of 



and 



'X' 



assuming the mean £a 



1 value for Segiie 
^, while Figure 
fi^fi^ (left) an 



1 in Eq. OL for the channel xX ~^ '- 
shows the same for the channels xX 
XX ^ 00 ^ /i+/^-/i+/z- (right). 

Figure [3] gives an indication of how well IceCube can 
probe the "best-case" scenario of dark matter annihila- 
tions directly into neutrinos. The right plot in this figure 
shows the expected neutrino exclusion limits from Ice- 
Cube observations of Segue 1 for this channel. 

Figure [4] shows that if dark matter with m^ >1 TeV 
explains the PAMELA and/or Fermi excesses, one could 
expect a few muon signal events from Segue 1 for the 
first of these channels, but perhaps only one for the 
annihilation channel through an intermediary 0. Note 
that the PAMELA/Fermi best fit regions assume a lo- 
cal dark matter density of 0.43 GeV/cm^ 23 — if the 
local dark matter density is lower, then a few more neu- 
trinos could be expected. Moreover, if the >Cann value 
is larger than the mean, then also a few more could be 
expected. The number of expected neutrinos could also 
be larger in the case of Sommerfeld enhanced models, 
where the low velocity dispersion of the dwarf galaxy may 
lead to a larger dark matter annihilation cross-section 
in the dwarf galaxy compared to locally (see below). 
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Figure [5] shows the expected neutrino exclusion hmits 
from IceCube observations of Segue 1 for the same chan- 
nels. For both annihilation channels, the ACT prospects 
are better than those from IceCube for low dark matter 
masses; for high dark matter masses (above the preferred 
PAMELA/Fermi regions) both are comparable. How- 
ever, the observing time for ACTs is of order a few days 
whereas for IceCube it is years. 

We note that neutrinos from dwarf galaxies are less 
sensitive probes of dark matter annihilation to charged 
leptons than are observations of neutrinos from the 
Galactic center region [5D] . However, observations of dif- 
ferent dwarf galaxies may be stacked, increasing the sen- 
sitivity by a factor of a few, and any resulting neutrino 
signal together with a corroborating gamma-ray signal 
would be very strong evidence for dark matter annihila- 
tion. 



Constraining Sommerfeld enhanced models 

In order for annihilating dark matter to explain the 
PAMELA and/or Fermi/HESS e~ /e^ excesses, the anni- 
hilation cross-section must be 0(100— 1000) times larger 
than the thermal WIMP freeze-out cross-section, which 
we take to be 



{crVrc\)i.o. 



2.5 X 10" 



"26 cm^s" 



(17) 



If the WIMP was produced non-thermally in the early 
Universe (for example, from the decay of a heavier par- 
ticle), it can naturally have a large annihilation cross- 
section [301 IHTVt63) . However, even if the WIMP was in 
thermal equilibrium in the early Universe, several mech- 
anisms allow the cross-section to be large today while be- 
ing much smaller at freeze-out P^HTOl [MH55J . The Som- 
merfeld enhancement [l4Hl6l [65H68] increases the cross- 
section at low dark matter velocities, v, so that the cross- 
section can be much larger today in the Milky- Way halo, 
where v ~ 10~^, than during freeze-out, where v ~ 0.3. 
Since the dark matter velocity in dwarfs is f ^ 10~^, 
there may be an additional order of magnitude enhance- 
ment in the annihilation cross-section in dwarfs (see also 
[SHHTT]). This additional enhancement is not guaran- 
teed, however, since the Sommerfeld effect may saturate 
for very low v. (In this paper, we do not consider any ad- 
ditional enhancement from the even colder substructure 
within the dwarf galaxy, but see for example [72j.) 

The Sommerfeld enhancement occurs when the dark 
matter particle x couples to a new mediator (/>, with mass 
TO0 ^ m^ , that produces an attractive potential between 
the dark matter particles. The 4> in turn couples to Stan- 
dard Model matter — in particular, we imagine that 
can decay to e^e~ or fJ-'^ fJ-~ , so that dark matter annihila- 
tion XX ^ 00 -> e+e"e+e" or XX^ H^ M^/""M^Ai", 
or a combination, can produce the cosmic-ray excesses (in 
many models, and depending on m^, (j) decays will also 
produce light hadrons). In this section, we show that 
gamma-ray observations of Segue 1 disfavor low masses 



for 4>, m^ ~ 0(0.1 GeV), assuming the mean value for 
£ann, although the uncertainty in £ann does not allow 
one to conclusively rule out such low masses. 

To calculate the Sommerfeld enhancement, we assume 
that the dark matter couples to a mediator with cou- 
pling strength A, and we consider two cases. The first 
case assumes that the dark matter interaction with is 
purely elastic^ while the second case assumes that the in- 
teraction is purely inelastic. The latter case occurs in the 
presence of an additional excited dark matter state x* so 
that the interaction with is only off-diagonal (i.e. x 
couples only inelastically to x*0j ^nd not elastically to 
x4>)- We denote the mass splitting between x and x* by 
5. The Sommerfeld enhancement in the presence of an 
excited state has been investigated in detail in ^75] , where 
it was found that generically it predicts larger enhance- 
ments compared to the elastic case. Excited states occur 
in many models, including the "inelastic dark-matter" 
[73] and "exciting dark matter" models |75j . 

For the elastic case, we model the interaction with an 
attractive Yukawa potential of the form (following [15] ) 



r 



(18) 



where ajj — X"^ /Att. The wave- function of two annihilat- 
ing dark matter particles in the non-relativistic limit is 
found by solving the Schrodinger equation, 

V'"(r) - m^V{r)jp{r) + v^m\%l:{r) = 0, (19) 

with the boundary condition ■(/''(oo) = ivmy.'ip{oo) and 
with ^(0) arbitrary. The total annihilation cross-section 
is then given by 



{<7V 



rcU 



So{aVrc\)o 
2 



ip{oo) 



^(0) 



(cWrcOo, 



(20) 



where v = Wroi/2 is the dark matter particle's velocity 
in the center-of-mass frame, Sq denotes the Sommerfeld 
enhancement in the elastic case, and 



(crWrcl)o 



"D 



(21) 



is the cross-section for XX ~^ 4>4> annihilation. This 
cross-section is modified if there are additional states into 
which X can annihilate, but we will not consider this pos- 
sibility (however see [751|77|). In the limit v <^ 1 and 
TO0 ^ TOj^, Eq. (19) may be solved exactly to give 



Sq 



a^TT 



V 



1 



-OLD-K JV 



(22) 



which clearly shows the 1/w dependence of the cross- 
section. Away from this limit, the solution to equation 
( [19| ) must be numerically integrated; it has resonances 
and saturates when v '^ m^/m^. Note that for ajj <^ v, 



Eq. ( 22 ) can be expanded to give 
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FIG. 

m^ = 0.1 GeV and dark matter mass splitting 5 = (upper left), m^ = 1 GeV and 5 = (upper right), m^ = 0.25 GeV 
and (5=1 MeV (lower left), and m^ = 1.0 GeV and (5 = 1 MeV (lower right) assuming the dark matter velocity dispersion of 
the Milky- Way halo (lower purple line) or Segue 1 (upper blue line). Here 5 = refers to dark matter without excited states, 
while 5 7^ refers to the splitting between the dark matter ground state and an excited state. The coupling an between the 
dark matter and (;A is a function of m^ and has been set to the value that gives the correct thermal relic abundance. Note that 
/3 = v/c. Limits from MAGIC observations of Segue 1 are shown in black solid or dashed lines, where the solid line uses the 
mean value of £ann and the upper (lower) dashed line uses the 2a lower (upper) bound on £ann for Segue 1 (see Eq. ([3|). Green 
and orange regions indicate the best-fit regions to explain the PAMELA and Fermi e^/e~ signals for (j> — >■ e^e" (upper plots) 
and for 4> — >■ fJ-'^ n~ (lower plots) and are adapted from [2D] after rescaling from a local dark matter density of 0.3 GeV/cm^ to 
0.43 GeV/cm'' [47]. Note that the purple lines intersect, or at least come close to, the PAMELA/Fermi preferred regions. This 



figure shows that a wide range of values for m^ can be probed and that, in particular, rritf, < 
current constraints (upper left plot). 



0.1 GeV is slightly disfavored by 



A useful analytical formula [68l [78] that closely 
matches the numerical solution can be obtained by us- 
ing the Hulthen potential V(r) = —aDm^e~^"^'>- /{I — 
e"'""*'*), which has similar behavior to the Yukawa po- 
tential for r — >■ and r — >■ cx). The solution is given 
by 
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We will use this analytic formula in the calculation of the 
Sommerfeld enhancement in the elastic case as described 
below. 

In order to calculate the Sommerfeld enhancement in 
the Milky- Way halo and in Segue 1 today, we proceed as 
follows (see also [76l Ull [79] ) . First we use Eqs.j2l]) and 
( 23 1 to set the annihilation cross-section in Eq. ( 20 ) equal 



(24) to the cross-section that gives the correct thermal relic 
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abundance during freeze-out, Eq. (17). This allows us to 



capture the Sommerfeld effect during freeze-out, where, 
although much less important than in the Milky- Way and 
dwarf halo now, it can still enhance the dark matter anni- 
hilation cross-section by an 0(1) amount [751 [771 [5nH52] . 
For a given m^, this allows us to calculate ao- The 
precise value of m^ is not very important, and the use 



of the simple formula Eq. ( 23 1 as opposed to its more 
complicated variations is sufficient. Moreover, we use 
V ^ 0.1 (as opposed to the canonical value v ^ 0.3) 
for the velocity of the dark matter during freeze-out and 
have checked that the resulting value of ao as a func- 
tion of m^ agrees to within about 10% with a detailed 
calculation of the thermal relic abundance in [76]. The 
resulting value of ajj may now be used to calculate the 
Sommerfeld enhancement in the Milky- Way and Segue 1 
halo. In particular, we take the dark matter distribution 
to be thermal so that the thermally averaged Sommerfeld 
enhanced annihilation cross-section is given by ^75] 
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(27) 



where N = erf(z/-\/2) — ^/2/n ze^^ /^ with z = Wcsc/^Oi 
fosc is the dark matter escape velocity, and Wmax is the 
maximum relative dark matter velocity in the laboratory 
frame, which wc will simply set equal to Wesc • Here we 



use the analytic formula for 5*0 given in Eq. (24), and 
take for the Milky- Way [53] 



vo =i 210 kni/s. 



525 km/s 



(28) 



(these values are close to other values in the literature, 
e.g. IMIIHS], and allow us to compare our results to those 
in "76]). For Segue 1, the average one-dimensional dark 
matter velocity dispersion within a 0.25° radius from the 
center is 4.5±1.8 km/s. This was obtained using Eq. (|4|. 
The relevant parameters we use for Segue 1 are 

vq ~ V 2 X 4.5 km/s ~ 6.4 km/s, Uesc ^ I'max ^ 10 km/s. 

(29) 
We note that small variations in the value of Wcsc and 
I'max for Segue 1 do not qualitatively impact the results. 
In Figure p] (upper two figures), we show the result- 
ing Sommerfeld enhanced annihilation cross-section as 
a function of m^ for m^ = 0.1 GeV (upper left) and 
m^ = 1 GeV (upper right) for dark matter annihilating 
in the Milky- Way halo (purple line) and in Segue 1 (blue 
line). We superimpose these on the PAMELA/Fermi 
preferred region adapted from [20], for the case that <j) 
decays exclusively to e+e^, after rescaling their values 
from a local dark matter density of 0.3 GeV/cm'^ to 0.43 
GeV/cm^ [12|M|- Note that the purple hue does inter- 
sect the PAMELA/Fermi preferred region, as it must if 
it is to explain this data. The gamma-ray constraints 



from MAGIC observations of Segue 1 are shown in black 
lines, with the solid line using the mean value of £ann, 
while the upper dashed shows the conservative limit with 
'Cann — 2(7, and the lower dashed line shows the optimistic 
limit with >Cann+2cr. We see that the additional Sommer- 
feld enhancement expected in Segue 1 (indicated by the 
blue line) allows us to constrain this model. In particu- 
lar, low values of m^ are disfavored assuming the mean 
value for £ann, although they cannot be ruled out due 
to the uncertainty in >Cann- Since the Sommerfeld en- 
hancement saturates when v < m^/m^, larger m^ are 
less constrained. However, even for larger m^ '^ 0(1 
GeV) a sizeable fraction of the parameter space lies near 
a resonance region, and is constrained. 

For the elastic case and for decaying exclusively to 
fi^ fi^ , the ao required by thermal freeze-out is too small 
to give a large enough Sommerfeld enhanced annihilation 
cross-section to explain the PAMELA/Fermi regions 76j . 
Inelastic interactions increase the cross-section and alle- 
viate the discrepancy [731 [77] . In Figure [6] (lower two 
plots), we show the resulting Sommerfeld enhanced an- 
nihilation cross-section as a function of m^, for 6 = 1 
MeV for m^ = 0.25 GeV (lower left) and m^ = 1 GeV 
(lower right) for dark matter annihilating in the Milky- 
Way halo (purple line) and in Segue 1 (blue hue). We 
review the details of the Sommerfeld enhancement in the 
presence of excited states in Appendix IA] (note that we 
ignore any complications in the freeze-out calculation of 
having additional excited states). For the case in which 
(j) decays exclusively to n'^n~, we superimpose these on 
the PAMELA/Fermi preferred region adapted from |20j . 
again rescaling their results to a local dark matter den- 
sity of 0.43 GcV/cm^ (note, however, that the Sommer- 
feld cross-sections are independent of the (j) decay modes) . 
The gamma-ray constraints from MAGIC observations of 
Segue 1 are again shown in black lines, and we find that 
for both small and larger values of m^, they constrain 
part of the relevant parameter space. 

The purple lines in the two lower plots of Figure [6] lie 
just below the Fermi region. This small discrepancy can 
be alleviated in many ways. For example, the local dark 
matter density may be larger than the value assumed 
here, 0.43 GeV/cm"^. Moreover, taking into account sub- 
structure which enhances the positron and electron flux 
can lower the required cross-section by another factor 
of two \E^, f57] . Uncertainties in the astrophysical back- 
grounds and in the cosmic-ray propagation can further 
lower the required cross-section and remove the discrep- 
ancy (see e.g. [HHllHS])- In addition, the PAMELA/Fermi 
regions assume that </> annihilates exclusively to fj,'^ ^~ . 
In many models, however, will decay also at least partly 
to e"'"e~, which lowers the required cross-section further, 
while still providing an excellent fit to the data. This is 
the case, for example, in which dark matter is charged un- 
der a new GeV-scale Abelian force that kinetically mixes 
with hypercharge, see e.g. [ISHISl 175 1 177] . 

The next generation of ACTs such as AGIS [lOOj 
and CTA '90j will improve the current prospects by 
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roughly an order of magnitude, and will significantly 
constrain Sommerfeld enhanced models that explain 
the Fermi/PAMELA signals. The current and future 
prospects from gamma-ray observations of dwarf galax- 
ies are complementary to those coming from current 
constraints of WMAP observations of the Cosmic Mi- 
crowave Background and prospective constraints from 
Planck 191-94:. 



IV. CONCLUSIONS 

The Milky Way satellite Segue 1 is the least luminous 
and most dark-matter dominated galaxy known. Us- 
ing a new sample of member stars, we have derived a 
new value for the line-of-sight integral of the dark mat- 
ter halo density-squared, >Cann- Its mean value is simi- 
lar to the same quantity for Draco and larger than that 
for all other known dwarfs. With this new value, we 
have derived robust upper limits on dark matter anni- 
hilation cross-sections from gamma-ray flux upper lim- 
its from Fermi and the Imaging Atmospheric Chcrenkov 
Telescope (ACT) MAGIC, and have shown that these 
constraints may in fact be much stronger given the sig- 
nificant uncertainty remaining in the Segue 1 dark matter 
distribution. For a mean value of the Segue 1 density- 
squared integral, the constraints we present are more 
stringent than those obtained from any other Milky Way 
dwarf spheroidals to date (see e.g. [5T1[^[5B] ). 

We showed that Fermi and ACTs are complementary 
in the dark matter models they can probe, with Fermi 
having an advantage over ACTs for dark matter channels 
that produce softer photons as well as for lower dark 
matter masses < 0(100 GeV). ACTs have an advantage 
over Fermi in constraining cross-sections for larger dark 
matter masses > 0(100 GcV). 

The flux limits do not rule out a dark matter explana- 
tion for the PAMELA/Fermi/HESS e+/e- signals, but 
probe the relevant regions. With its low dark matter ve- 
locity dispersion. Segue 1 also serves as a useful labora- 
tory for constraining Sommerfeld enhanced models that 
have been invoked to explain these signals. We find that 
current constraints from MAGIC slightly disfavor light 
mediators of mass 0(0.1 GeV) if we take the mean value 
for £ann- Howcver, due to the uncertainty in the dark 
matter density in Segue 1, such light force carriers can- 
not at this stage be ruled out. 

In addition, we show that under favorable circum- 
stances neutrino observatories such as IceCube also have 
the potential to detect a dark matter signal from Segue 
1, providing an important cross-check, for some channels, 
if a signal is detected in the gamma-rays. 

Since Segue 1 lies in the Northern Hemisphere, it is 
an excellent target for the VERITAS and MAGIC, and 
for IceCube. Observations by MAGIC (and future obser- 
vations by VERITAS) are able to set interesting and ro- 
bust constraints on dark matter annihilation, with rather 
modest observing times directed at Segue 1. Moreover, 
the ACT constraints presented here may be considerably 



improved if these instruments are able to lower their en- 
ergy thresholds, observe over longer times, and perform 
an analysis specifically aimed at models that explain the 
e+/e^ signals. 

Though our analysis has focused on the galaxy Segue 
1, prospects are promising for strengthening the limits 
we have presented. Given selection effects and biases in 
current surveys for MW satellites, it is probable that the 
population of nearby, dark-matter dominated satellites 
will increase in the near future |97j . Any signal from 
one dwarf can thus be corroborated by looking at oth- 
ers. Moreover, the significance of a signal or of a limit 
can be improved by stacking together the observations of 
several dwarf galaxies. In addition, the next generation 
of ACTs such as AGIS and CTA will improve the cur- 
rent prospects by roughly an order of magnitude. These 
improvements will allow a definitive probe of many dark 
matter annihilation models. 
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Appendix A: Sommerfeld enhancement in the 
presence of an excited dark matter state 

In this brief Appendix, we summarize the formulas 
needed to calculate the Sommerfeld enhancement in the 
presence of an excited state [73]. The Sommerfeld ef- 
fect can be calculated numerically, however |73| presented 
a computationally less intensive semi- analytical formula, 
which we review here. We refer the reader to [73| for all 
the details, and mostly use the same notation. 

We take the mass-splitting between the dark matter 
ground state and excited state to be S, and define 



V m0 
, £0 — : e<5 




(Al) 



For El,, £0, e^- ^ 1 and for 5 < max(a£)7TT-0, TOx(w/c)^), the 
Sommerfeld enhanced annihilation of the ground-state is 
given by 
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sh(^)-cos(^7Jf3l2+2e_) 

cosh U(e„ + y=I|+;!")j-cos(2e_) 



Eu < £5, 



eu > £5- 



(A2) 



Here 




- 2a*zJ - 4zi - J ^' dr Vl-A-I 



(A3) 
(A4) 



where r^/ and zm are defined by 



rM 



zm — max 



and A_ is given by 



A-M 



{pi) 

( 4 4 ] 

V6V4' 16^4;' 



1 „\2 e-2^<*'- 

2 -^ 



(A5) 
(A6) 

(A7) 



r 



In the lower two figures in Figure [6j we show the Som- 
merfeld enhanced annihilation cross-section for 5 = 1 
MeV. Note that here we did not thermally average the 
cross-section, which is computationally very intensive, 
but rather took 150 km/s for the Milky- Way dark matter 



velocity and 4.5 km/s for that of Segue 1. For the region 
of interest around m^ ~ 0(1 TeV), we checked that the 
resulting cross-section is ~ 20% lower than the thermally 
averaged cross-section. 



[1] V. Belokurov et al. (SDSS), Astrophys. J. 654 

(2007), astro-ph/0608448. 
[2] M. Geha, B. Willman, J. D. Simon, L. E. Strigari, 

Kirby, D. R. Law, and J. Strader, Astropliys. J. 

1464 (2009), 0809.2781. 
[3] M. Niederste-Ostholt et al. (2009), 0906.3669. 
[4] J. D. Simon et al. (2010), 1007.4198. 
[5] M. Geha, TeV Particle Astrophysics, 2009. 
[6] L. Bergstrom, T. Bringmann, M. Eriksson, 

M. Gustafsson, Phys. Rev. Lett. 94, 131301 (2 

astro-ph/0410359. 
[7] J. F. Beacom, N. F. Bell, and G. Bertone, Phys. 

Lett. 94, 171301 (2005), astro-ph/0409403. 
[8] A. Birkedal et al. (2005), arXiv:hep-ph/0507194. 
[9] O. Adrian! et al. (2008), 0810.4995. 
[10] A. A. Abdo et al. (The Fermi LAT) (2009), 0905.1 
[11] F. Aharonian et al. (H.E.S.S.), Phys. Rev. Lett. 

261104 (2008), 0811.3894. 
[12] H. E. S. S. C. F. Aharonian (2009), 0905.0105. 



, 897 

E.N. 
692, 


[13] 
[14] 
[15] 
[16] 




[17] 


and 


[18] 


!005), 


[19] 


Rev. 


[20] 




[21] 


0025. 
101, 


[22] 




[23] 
[24] 



0. Adrian! et al. (2008), 0810.4994. 
M. Cirell! et al. (2008), 0809.2409. 

N. Arkani-Hamed et al. (2008), 0810.0713. 

M. Pospelov and A. Ritz, Phys. Lett. B671, 391 (2009), 

0810.1502. 

1. Cholis, D. P. Finkbeiner, L. Goodenough, and 
N. Weiner (2008), 0810.5344. 

L Cholis, G. Dobler, D. P. Finkbeiner, L. Goodenough, 
and N. Weiner (2008), 0811.3641. 

L. Bergstrom, J. Edsjo, and G. Zaharijas (2009), 
0905.0333. 

P. Meade, M. Papucci, A. Strumia, and T. Volansky 
(2009), 0905.0480. 

G. D. Martinez, J. S. Bullock, M. Kaplinghat, L. E. Stri- 
gari, and R. Trotta, JCAP 0906, 014 (2009), 0902.4715. 
R. Essig, N. Sehgal, and L. E. Strigari, Phys. Rev. D80, 
023506 (2009), 0902.4750. 

P. Scott et al., JCAP 1001, 031 (2010), 0909.3300. 
M. Perelstein and B. Shakya (2010), 1007.0018. 



14 



[25] Q. E. Minor, G. Martinez, J. Bullock, M. Kaplinghat, 

and R. Trainor (2010), 1001.1160. [59 

[26] L. E. Strigari et al., Astrophys. J. 678, 614 (2008), 

0709.1510. [60 

[27] J. F. Navarro et al. (2008), 0810.1522. 

[28] M. Regis and R UUio (2009), 0907.5093. [61 

[29] P. Colin, V. Avila-Reese, and O. Valenzuela, Astrophys. 

J. 542, 622 (2000), astro-pli/0004115. [62 

[30] L. M. Widrow and J. Dubinski, Astrophys. J. 631, 838 

(2005), arXiv:astro-ph/0506177. [63 

[31] A. A. Abdo et al., Astrophys. J. 712, 147 (2010), 

1001.4531. [64 

[32] B. Willman et al. (2010), 1007.3499. 
[33] R. A. Ibata, G. F. Gilmore, and M. J. Irwin, Mon. Not. [65 

Roy. Astron. Soc. 277, 781 (1995), astro-ph/9506071. [66^ 

[34] G. D. Mack et al., Phys. Rev. D78, 063542 (2008), 

0803.0157. [67 

[35] N. F. Bell and T. D. Jacques (2008), 0811.0821. 
[36] J.-F. Fortin, J. Shelton, S. Thomas, and Y. Zhao (2009), [ 

0908.2258. 
[37] N. Fornengo, L. Fieri, and S. Scopel, Phys. Rev. D70, [69 

103529 (2004), hep-ph/0407342. [70^ 

[38] Y. Kuno and Y. Okada, Rev. Mod. Phys. 73, 151 (2001), [71 

hep-ph/9909265. [72 

[39] J. Mardon et al. (2009), 0901.2926. [73 

[40] G. Kane, R. Lu, and S. Watson (2009), 0906.4765. [74 

[41] P. Gondolo et al., JCAP 0407, 008 (2004), astro- 

ph/0406204. [75 

[42] http://www.physto.se/~edsjo/darksusy. 
[43] V. Barger, W.-Y. Kcung, G. Shaughnessy, and A. Tre- [76 

gre, Phys. Rev. D76, 095008 (2007), 0708.1325. 
[44] P. Sandick, D. Spolyar, M. R. Buckley, K. Freese, and [77 

D. Hooper, Phys. Rev. D81, 083506 (2010), 0912.0513. 
[45] C. Amsler et al. (Particle Data Group), Phys. Lett. [78 

B667, 1 (2008). 
[46] E. Rcsconi (IceCube Collaboration), Nucl. lustrum. [79 

Meth. A602, 7 (2009), 0807.3891. 
[47] P. Salucci, F. Nesti, G. Gentile, and C. F. Martins [80 

(2010), 1003.3101. [81 

[48] V. Barger and R. Phillips, Collider Physics (Addison- 

Wesley (Redwood City, USA), 1987). [82; 

[49] K. Winter, Neutrino Physics (Cambridge University 

Press, 2000). [83 

[50] J. Ahrens, J. N. Bahcall, X. Bai, R. C. Bay, T. Becka, 

K.-H. Becker, D. Berley, E. Bernardini, D. Bertrand, [ 

D. Z. Besson, et al., Astroparticle Physics 20, 507 

(2004), arXiv:astro-ph/0305196. [85 

[51] M. C. Gonzalez-Garcia, F. Halzen, and M. Maltoni, 

Phys. Rev. D71, 093010 (2005), hep-ph/0502223. [86 

[52] M. Honda, T. Kajita, K. Kasahara, S. Midorikawa, 

and T. Sanuki, Phys. Rev. D75, 043006 (2007), astro- [87 

ph/0611418. 
[53] J. Ahrens et al. (IceCube), Astropart. Phys. 20, 507 [ 

(2004), astro-ph/0305196. 
[54] P. Wang (Fcrmi-LAT), Cosmology in Northern Califor- [ 

nia Meeting, 2009. [90| 

[55] C. Farnier (Fermi-LAT), Roma International Confer- 
ence on Astro-Particle Physics, 2009. [91 
[56] S. Murgia (Fermi-LAT), TeV Particle Astrophysics, 

2009. [92 

[57] T. Jeltema (Fermi-LAT), TeV Particle Astrophysics, 

2009. [93 

[58] D.N. Castafio, Dark Matter 2010, UCLA, 

http://www.physics.ucla.edu/hep/dmlO/talks/ [94 



nietocastano.pdf. 

A. A. Abdo et al., Phys. Rev. Lett. 104, 091302 (2010), 

1001.4836. 

D. Spolyar, M. R. Buckley, K. Freese, D. Hooper, and 

H. Murayama (2009), 0905.4764. 

M. Kamionkowski and M. S. Turner, Phys. Rev. D 42, 

3310 (1990). 

T. Moroi and L. Randall, Nucl. Phys. B570, 455 (2000), 

hep-ph/9906527. 

G. F. Giudice, E. W. Kolb, and A. Riotto, Phys. Rev. 

D64, 023508 (2001), hep-ph/0005123. 

M. Ibe, H. Murayama, and T. T. Yanagida, Phys. Rev. 

D79, 095009 (2009), 0812.0072. 

A. Sommerfeld, Annalen der Physik 403, 257 (1931). 

J. Hisano, S. Matsumoto, and M. M. Nojiri, Phys. Rev. 

Lett. 92, 031303 (2004), hep-ph/0307216. 

J. March-Russell, S. M. West, D. Cumberbatch, and 

D. Hooper, JHEP 07, 058 (2008), 0801.3440. 

J. D. March-Russell and S. M. West, Phys. Lett. B676, 

133 (2009), 0812.0559. 

B. Robertson and A. Zentner (2009), 0902.0362. 
M. Kuhlen (2009), 0906.1822. 

M. Kuhlen, P. Madau, and J. Silk (2009), 0907.0005. 
J. Bovy, Phys. Rev. D79, 083539 (2009), 0903.0413. 
T. R. Slatyer, JCAP 1002, 028 (2010), 0910.5713. 
D. Tucker-Smith and N. Weiner, Phys. Rev. D64, 
043502 (2001), hep-ph/0101138. 

D. P. Finkbeiner and N. Weiner, Phys. Rev. D76, 
083519 (2007), astro-ph/0702587. 

J. L. Feng, M. Kaplinghat, and H.-B. Yu (2010), 
1005.4678. 

D. P. Finkbeiner, L. Goodenough, T. Slatyer, and 
N. Weiner, to appear (2010). 

C.-Y. Chen, F.-L. Lu, and D.-S. Sun, Physica Scripta 
76, 428 (2007). 

J. L. Feng, M. Kaphnghat, and H.-B. Yu, Phys. Rev. 
Lett. 104, 151301 (2010), 0911.0422. 
Q. Yuan et al., JCAP 0912, Oil (2009), 0905.2736. 
J. B. Dent, S. Dutta, and R. J. Scherrer, Phys. Lett. 
B687, 275 (2010), 0909.4128. 

J. Zavala, M. Vogelsberger, and S. D. M. White, Phys. 
Rev. D81, 083502 (2010), 0910.5221. 
X. X. Xue et al. (SDSS), Astrophys. J. 684, 1143 (2008), 
0801.1232. 

M. J. Reid, K. M. Menten, A. Brunthaler, and G. A. 
Moellenbrock (2009), 0902.3928. 

M. C. Smith et al., Mon. Not. Roy. Astron. Soc. 379, 
755 (2007), astro-ph/0611671. 

M. Kuhlen and D. Malyshev, Phys. Rev. D79, 123517 
(2009), 0904.3378. 

J. M. Cline, A. C. Vincent, and W. Xue, Phys. Rev. 
D81, 083512 (2010), 1001.5399. 

L. Stawarz, V. Petrosian, and R. D. Blandford, Astro- 
phys. J. 710, 236 (2010), 0908.1094. 
B. Katz, K. Blum, and E. Waxman (2009), 0907.1686. 
T. Bringmann, M. Doro, and M. Fornasa, JCAP 0901, 
016 (2009), 0809.2269. 

N. Padmanabhan and D. P. Finkbeiner, Phys. Rev. 
D72, 023508 (2005), astro-ph/0503486. 
S. Galli, F. locco, G. Bertone, and A. Melchiorri, Phys. 
Rev. D80, 023505 (2009), 0905.0003. 
T. R. Slatyer, N. Padmanabhan, and D. P. Finkbeiner 
(2009), 0906.1197. 
S. Galh et al. (2010), 1005.3808. 



15 



[95] J. Albert et al. (MAGIC), Astrophys. J. 679, 428 

(2008), 0711.2574. 
[96] V. Acciari et al. (The VERITAS) (2010), 1006.5955. 
[97] B. Willman (2009), 0907.4758. 
[98] We note that including the two RR Lyrae variable stars 

and two probable binaries has minimal effect on our 

results. 
[99] We note that the cosmic-ray spectrum obtained from 

(j)'s decaying exclusively to e"'"e~ gives a slightly worse 



[100] 



fit to the PAMELA/Fermi spectra than if one, for ex- 
ample, assumes additional hidden sector showering, or 
4> decays to electrons and muons, or only to muons [5D]. 
We refer the reader to [77] for a detailed discussion and 
benchmark points consistent with the PAMELA/Fermi 
spectra and preferred regions in the cross-section and 
mass plane. 
www.agis-observatory.org 



